OTIS  FILE  COEy  AD-A147  855 


' 

t 

OFFICE  OF  NAVAL  RESEARCH 
Contract  N00014-80-C-0472 
Task  No.  NR  056-749 
TECHNICAL  REPORT  No.  57 


Theoretical  Study  of  Laser-Induced 
Surface  Excitations  on  a  Grating 

by 

Ki-Tung  Lee  and  Thomas  F.  George 


Prepared  for  Publication 
in 

Physical  Review  B 


Department  of  Chemistry 
University  of  Rochester 
Rochester,  New  York  14627 


November  1984 


Reproduction  in  whole  or  in  part  is  permitted  for  any 
purpose  of  the  United  States  Government. 


This  document  has  been  approved  for  public  release 
and  sale;  its  distribution  is  unlimited. 


4.  PERFORMING  ORGANIZATION  REPORT  NUMBER  IS) 


B.  MONITORING  ORGANIZATION  REPORT  NUMBER(f) 


UR0CHESTER/DC/84/TR-57 


(4  NAME  OF  PERFORMING  ORGANIZATION  IMl  OFFICE  SYMBOL 

Department  of  Chemistry  I  »t* *#*•*•> 
University  of  Rochester 


Sc.  ADDRESS  (City.  Slat*  and  ZIP  Cadi) 

River  Station 
Rochester,  New  York  14627 


Bp.  name  OF  FUNOING/SPONSORING 

8f?ice Ao^aval  Research 


Be.  AOORESS  (City.  Slate  and  ZIP  Coda) 

Chemistry  Program  Code  472 
Arlington,  Virginia  22217 


11  TT^Theoretical  Study  of  Laser-Induced  Surface  Excitations  on  a  Grating 


12.  PERSONAL  AUTMORIS)  T _ , _ _ .  Tl _ r  - _ 


.OFFICE  SYMBOL 
Ilf  applicable) 


7a.  NAME  OF  MONITORING  ORGANIZATION 

Office  of  Naval  Research 


7b.  AOORESS  fCMy.  Slate  and  ZIP  Coda) 

Chemistry  Program  Code  472 
Arlington,  Virginia  22217 


B.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 

Contract  N00014-80-C-0472 


1 10.  SOURCE  OF  FUNDING  NOS. 

PROGRAM 
ELEMENT  NO. 

NR  056-749 

PROJECT 

NO. 

TASK 

NO. 

Ki-Tung  Lee  and  Thomas  F.  George 


13a.  TYPE  OF  REPORT  . 

Interim  Technical 


13b.  TIME  COVERED 
FROM _  TO _ 


IB.  SUPPLEMENTARY  NOTATION 

Prepared  for  publication  in  the  Physical  Review  B 


18.  PAGE  COUNT 


IE  SUBJECT  TERMS  IContfnuc  on  raiwnw  Ifntcemm ry  and  Identify  By  Block  numBcr) 

GRATINGS  MAXWELL'S  EQUATIONS 

SyK£CS:!LX.SJIATI0NS  ARBITRARY  GRATING  SHAPES 

LASER- INDUCED  NEW  RESONANCE  PHENOMENON 


17. 


COSATI  CODES 


THEORETICAL  STUDY  OF  LASER-INDUCED  SURFACE  EXCITATIONS  ON  A  GRATING 


XI -Tung  Lee  end  Thomas  F.  George 

Department  of  Chemistry 
University  of  Rochester 
Rochester.  New  York  14627 


Laser-Induced  surface  excitations  on  a  grating  are  studied  In  terms  of  the 
solutions  to  Maxwell's  equations.  A  rigorous  theory,  derived  originally  for  a 
lamellar  grating.  Is  used  to  study  the  resonance  phenomenon  for  deep  gratings. 

The  generalization  of  the  square-well  grating  to  gratings  of  arbitrary  shapes  Is 
examined  numerically.  A  new  diffraction  anomaly  Is  seen  to  occur  when  the  grating 
depth  is  approximately  equal  to  half  of  the  wavelength  of  the  Incident  laser  radi¬ 
ation. 


I.  Introduction 

It  has  been  well  established  that  impinging  light  on  a  rough  metallic  surface 
may  lead  to  resonant  excitation  of  surface  plasma  oscillations,  which  play  an 

important  role  in  interesting  surface  phenomena,  such  as  surface-enhanced  Raman 

12  3 

scattering  *  and  laser-induced  periodic  pattern  deposition.  In  particular, 

the  spatial  oscillation  of  photochemlcally  deposited  metal  films  is  identified 
as  the  fingerprint  of  the  oscillation  of  a  surface  plasma  wave  (SPU).  These 
surface  waves  are  related  or  similar  to  those  found  in  Wood's  anomalies4,5  or 
Brewster's  waves,5  which  all  satisfy  the  same  mathematical  equation  but  with  media 
of  different  electric  properties.  Moreover,  all  these  waves  have  the  same  charac¬ 
teristics  in  that  they  propagate  along  the  surface  of  a  dielectric  medium,  with 
amplitudes  decaying  exponentially  with  increasing  distance  from  the  surface  into 
the  dielectrics  and  into  the  vacuum  in  contact  with  it.  In  the  second  quantiza¬ 
tion  terminology,  these  surface  waves  are  labelled  as  surface  polaritons.  Resonant 
excitation  of  these  quasi  stationary  modes  can  lead  to  orders  of  magnitude  of  en¬ 
hancement  of  molecular  processes  occurring  near  or  on  the  surface.  (It  has  been 

7  8 

suggested  by  Tsang  and  Klrtley  and  later  by  Mill  and  Weber  that  the  maximum 
value  of  electric  fields  near  gratings  may  be  limited  by  grating-induced  radiative 
damping.)  Thus,  a  detailed  understanding  of  the  formation  and  the  properties  of 
these  surface  waves  would  be  helpful  In  better  understanding  and  controlling  surface 
molecular  rate  processes  and  diagnosis  of  the  properties  of  surfaces.  (A  compre¬ 
hensive  review  of  the  diagnosis  of  surface  properties  can  be  found  in  Reference  9.) 

A  general  understanding  of  the  physical  properties  of  these  surface  waves  has 
been  reviewed  by  Fano,5  and  recent  advances  in  research,  both  theoretically  and 
experimentally,  have  been  reviewed  by  Agranovich  and  Mllls.^  These  quasi  stationary 
modes  are  controlled  by  the  following  elements: 


a)  the  wavelength  of  the  Incident  radiation. 


b)  the  angle  of  Incidence, 

c)  the  geometry  of  the  surface  roughness, 

d)  the  electric  properties  of  the  surface, 

e)  the  polarization. 

One  approach  to  the  theory  of  surface  plasma  oscillations  has  been  in  terms 

of  the  collective  motion  of  electrons  in  the  solld.^  It  has  been  found,  however, 

that  the  problem  of  resonant  excitation  of  SPW  can  be  treated  more  easily  as  an 

optical  problem  than  as  the  collective  motion  of  the  electrons.  In  fact,  both  the 
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problems  of  surface-enhanced  Raman  scattering  and  laser-induced  pattern  deposition 

13 

have  been  examined  by  Rayleigh's  perturbative  diffraction  theory. 

The  formation  of  these  evanescent  waves  originates  a  momentum  transfer  from 
the  grating  (roughness)  to  the  Impinging  waves.  Since  the  frequency  of  the  Impinging 
laser  radiation  Is  unaffected  by  the  roughness,  the  modulus  of  the  momentum  of  any 
wave  must  equal  the  one  in  the  vacuum.  Thus,  momentum  transfer  by  the  roughness 
can  only  change  the  direction  of  momentum.  This  means  that  when  the  tangential 
component  of  the  momentum  of  the  diffractive  wave. 


m 


in 


where  Is  the  tangential  component  of  the  Incident  wave,  d  is  the  spacing  of  the 
grating  and  m  Is  an  Integer,  Is  larger  than  the  modulus  of  the  Incident  momentum, 
k0,  then  the  normal  component  of  the  diffractive  wave. 


k 


n 


(Ko-k 


2  >1/2 
mt; 


t 


becomes  Imaginary.  This  corresponds  to  evanescent  waves  travelling  with  momentum 
kmt  along  the  surface  and  exponentially  damped  In  the  normal  direction.  Moreover, 
If  the  geometry  Is  chosen  so  certain  klnematlcal  conditions  are  met,  the  Incident 
photon  will  resonantly  couple  to  these  surface  polarltons. 


yi. .-*„**! it'*.'*  ',£•  ■•**■■■*.  «  i..cr8.r <-<*•,  t'ttV'.Art ♦*«•»*>  i>ti.  A  *.!*,• 
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The  effect  of  the  grating  depth  on  the  resonant  conditions  has  been  studied 

recently  by  Garcia^  and  Nevlere  and  Relnlsch.^  While  Garcia  found  that  the 

Intensities  of  the  reflectivity  and  the  diffraction  beams  have  their  minimum  at 

resonance,  but  that  the  reflectivity  becomes  a  maximum  for  a  particular  value  of 

is 

the  grating  depth,  Nevlere  and  Relnlsch  pointed  out  that  there  exists  an  optimum 
value  of  the  grating  depth  for  which  the  electromagnetic  (EM)  resonance  contribution 
to  the  enhancement  of  the  nonlinear  optical  process  Is  the  greatest,  and  the  optimi¬ 
zation  is  achieved  with  very  shallow  modulation.  More  recently.  Glass,  Weber  and 
Mills, utilizing  the  extinction  theorem  of  Tolgo,  Marvin,  Hill  and  Celli,^ 
studied  the  grating-induced  radiative  damping  of  the  surface  polariton  for  grating 
profiles  with  various  shapes  and  depths,  and  compared  the  results  with  those 
of  perturbation  theory.  Good  agreement  was  found  for  those  systems  investigated. 

All  of  the  calculations  mentioned  above  fall  in  the  regime  that  the  grating  ampli¬ 
tude  is  much  smaller  than  the  Incident  wavelength.  In  this  paper,  we  employ  the 

ID 

square-well  grating  theory  of  Sheng,  Stepleman  and  Sandra,  °  which  Is  based  on  the 

19 

formalism  of  a  stratified  medium,  to  study  a  new  kind  of  resonance  condition  for 

a  grating  depth  comparable  to  the  Incident  wavelength. 

20 

Recently,  various  numerical  techniques  have  been  developed  for  the  solution 
of  diffraction  problems.  Although  most  of  these  are  rigorous  in  formulation,  they 
exhibit  various  degrees  of  effectiveness  In  actual  calculations.  Very  often 

21 

numerical  difficulties  are  found  either  when  the  grating  depth  becomes  too  large 
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or  when  the  conductivity  of  the  grating  Is  large.  The  square-well  grating  approach 
has  been  demonstrated  to  be  numerically  stable  even  In  the  regime  of  good  conductivity 
and  also  to  work  well  for  the  case  of  a  deep  grating.  However,  the  method  Is 

restricted  to  a  lamellar  grating.  We  have  suggested  a  generalization  of  the 

22 

square-well  grating  theory  to  gratings  of  arbitrary  shape.  In  this  paper,  the 
numerical  applicability  of  this  generalized  theory  Is  examined.  The  diffraction 


from  gratings  of  arbitrary  shape  has  attracted  a  lot  of  attention  recently. 

23 

Moharan  and  Gaylord  have  applied  their  coupled-wave  theory  to  a  multilayered 

25 

grating  formalism,  while  Suratteau,  Cadllhac  and  Petit  formulated  a 
"Multi-Step  Lamellar  Grating"  (MSLG)  method  by  extending  the  work  of 
Gotten  and  coworkers.  The  MSLG  method,  which  has  only  been  applied  to 
lossless  dielectrics.  Is  believed  to  be  closely  related  to  the  present  work. 
In  Section  II,  the  theory  of  the  square-well  grating  and  Its  generalization 
to  a  multilayered  grating  Is  briefly  reviewed.  The  results  and  discussion 
are  given  In  Section  III,  and  a  summary  Is  presented  In  Section  IV. 
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II.  Theory 

A.  Square-Well  Grating 


The  square-well  grating  formalism  Is  based  on  a  stratified  medium  model. 


17 


and  a  detailed  account  of  the  derivation  can  be  found  in  a  paper  by  Sheng  et  al. 
Our  main  purpose  in  this  part  of  Section  II  is  to  set  up  the  notation  for  the 
generalization  of  the  square-well  grating  to  a  multilayered  grating. 

Following  Rayleigh's  approach,  the  entire  space  Is  separated  Into  three 
regions:  the  vacuum  region  (Region  I),  the  dielectric  region  (Region  III)  and  the 
corrugation  region  (Region  II).  The  geometry  and  the  coordinate  system  of  the 
lamellar  grating  are  shown  In  Fig.  1.  In  Regions  I  and  III,  the  EM  fields  are 
expressed  In  terms  of  Rayleigh  waves. 
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y1  «  exp[1kQ(sineix  -  cos6.,z)]  "  I  Rn  *xP^ko^Ynx+  ^’T2)^2*]) 


(2) 


and 


f111  »  l  T„  exp{1k0[Yn*  -  (c-y‘),,cz]). 


2\l/2, 


(3) 


For  p-wave  scattering,  y  Is  the  y-component  of  the  magnetic  vector,  kQ  *  2rr/A 
Is  the  free-space  wave  vector  of  the  Incident  laser  radiation,  e  Is  the  dielectric 
constant  of  the  grating,  Rp  and  Tn  are  the  amplitudes  of  the  n-th  reflected  and 
transmitted  diffracted  orders,  0^  Is  the  angle  of  Incidence,  and 


Yn  *  sin  0^  ♦  nA/d, 


(4) 
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where  d  Is  the  period  of  the  grating.  (The  square  root  with  positive  Imaginary 
part  will  always  be  chosen  In  this  paper.)  The  field  In  Region  II  satisfies16 


»!4i  ♦  ^4L  +  e(x)k2*” .  asii  ^  i«[e(K)]. 


where  e(x)  «  c  for  |x-  nd|  £rd/2,  and  1  otherwise  (r  Is  a  number  between  0  and  1). 


At  this  stage,  the  problem  becomes  Identical  to  that  of  a  periodically  stratified  medi 


with  a  piecewise  constant  e.  Thus,  the  field  In  Region  II  can  be  expressed  In  gener 


*  *  l  X.(x)[A.exp(1ABz)  +  B.exp(»1A.z)], 


where  A*  are  the  eigenvalues  satisfying  the  transcendental  equation 


l  Q.  0. 

co$(lc0d  sine.,- )  +  cos  (B^rdjcosfa^  l-r)d]  ♦  £{eip+— ) 


£  £ 


X  sln(g^rd)  s1n[o4(l-r)d]  •  0, 


and  the  associated  eigenfunctions.  X.(x)  are  given  by 


hM 


cos[eJl(x+  -^j]  +  ivj  ^2  sinCe^x*  *£)],  |*|  <  ^ 


Ul  cos[at(x-^)]  ♦  iv'^sInCajlx-^)],  ^<x<(l-£)d 


where  X.(-rd/2)  Is  normalized  to  1,  with 


Vj  *  [exp(ikQd  sine^  -  M^/N^, 


_  1  H 


V1  *  ^-7  s1n(0trd)  ♦  Vj  cosU^rd), 


U 1  «  cosU^rd)  +  ivj  ~  slnte^rd) 

£ 


Bff 

=  cosjB^rd)  cosfo^O-rJd]  -e“ •  sin[aA(l-r)d]  slnfB^rd),  (12) 

£ 

Hl  =  ik0^a  cos<^rd)  s1n[aA(l-r)d]  +  sinfB^rd)  costc^O-rJd],  (13) 

A  £> 

“t  •  <ko  -  4>1/2  <14> 

Bt  =  (tk2  -  A2)1/2.  (15) 

The  remaining  task  is  to  determine  the  expansion  coefficients,  Rp,  Tn, 

A^  and  B^,  by  matching  the  analytic  solutions  of  Maxwell's  equations  at  the 
boundaries;  namely,  the  continuity  of  tangential  E  and  H  fields  at  z=0  and 
z  =  -h.  This  gives  four  simultaneous  equations  which  are  valid  for  all  x: 

l  exp(ik0Ynxj)(6n0-Rn)  =  IX^XjKA^+B*), 

n  x, 

Ixa(Xj)[AJlexp(-iA]lh)  +  B£exp(i  A^h)]  « 

I  exp{ik0[Ynx.+  (e-y2)1/2h]}Tn, 
n 

I  -  fk0[cos6,6n0  ♦  (l-Y2),/2R„e*p(ikoVj)]  ■ 
n 

l '  4  lxj>W  B»J 

„  MXi) 

l  iA^expt-IA^h)  -  B£exp(iA£h)]  * 

9* 

FI  -'ko<e-f»>1/2Tne*',{1koVj*  U-T*»,/!tH}- 


(16) 


(17) 


08) 


09) 


Employing  the  point  matching  method,  we  can  rewrite  Eqs.  (16)  -  (19)  In  matrix 
notation  as 


(20) 
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1 

r"1  x(r  a  +  r_1B)  *  at, 

ss  m  -  •  ** 

(21) 

-n(R+  o)  =  r“1n(A-  b), 

;  '  *  a  ~ 

(22) 

1 

r"1  a(ZA-Z-1B)  “-5  AT, 

«  *  a~  z  -  a  a- 

(23) 

g||  where 

Pjn  *  exP(ik0Ynxj)» 

(24) 

i 

XjA  “  Xsjxj^’ 

(25) 

1 

h 

fij£,  =  e  Vxj)* 

(26) 

1 

lH  =  exp(-iA^h)6j^, 

(27) 

1 

nnj  *  k0(,"'n>,/V 

(28) 

i 

4nj  *  «PClk0ve-Y^),/2h:«„j, 

(29) 

E 

£  .  .  I!®  (e.v2)V2{ 

Sj  e  'e  V  °nj* 

(30) 

D  is  defined  as  Dm  *  6,^.  After  some  manipulation*  we  have  the  expressions  for  R  and 
T  in  matrix  form  as 


with 


5  ■  (®-J)‘I(®  +  J)D 

s  *  -{[r’nter'Mr  +  (r',n-5r'1*)E',]n',r)'1 

*  »  *  SS  a  S  Z  X  ZZ  Z  Z  a  a 

x  [(r_1n  ♦  $r"1x)r  -(r“1n-er“1x)z“1]x’1r 
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(31) 


(32) 
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where 


k  ■  E“1(x‘1r-8’1r0  ♦  z(x'1r+n”1rc) - 


(34) 


For  h=0,  both  Z  and  A  reduce  to  identity  matrices.  In  this  case,  0  becomes 

*  *  ss 

-  5,  which  directly  implies 


-1 


R  =  (c+n)"  (c-n)D 


(35) 


and 


T  ■  D  -  R,  (36) 

which  are  the  correct  expressions  for  flat  surfaces.  In  the  case  of  a  very  deep 
grating,  h-*»,  and  for  with  an  imaginary  part  of  any  size,  then  0  can 

be  reduced  to  a  simple  form,  -flX-1,  and  the  partial  field  AT  goes  to  zero.  The 

s'* 

resonant  excitation  conditions  are  given  mathematically  by  the  zero's  of  the  expres¬ 
sion  (0HI)  (e-n)  and  physically  are  governed  by  the  five  parameters  given  in  Sec¬ 
tion  I.  Although  these  parameters  have  different  degrees  of  importance  in  the 
control  of  the  resonant  excitation  of  the  surface  plasma  oscillations,  one  can, 
with  appropriate  choice  of  these  five  parameters,  engineer  some  interesting  sur¬ 
face  phenomena.  Here,  we  focus  on  how  the  depth  of  the  square-well  grating  may 
introduce  such  phenomena. 


B.  Multilayered  Grating 

We  now  extend  the  square-well  grating  formalism  to  a  multilayered  grating 
theory  which  can  model  gratings  of  arbitrary  shape.  This  involves  the  extension 
of  the  concept  of  a  periodic  stratified  medium  to  a  two-dimensional  scheme.  A 
schematic  diagram  of  the  multilayered  grating  is  illustrated  in  Fig.  2,  where 


Region  II  is  divided  into  subregions.  In  each  subregion,  the  electromagnetic 
field  is  expressed  as  in  Eq.(6), 

-  Ex”(x)[A”exp(1^2) ♦  Bjexp(-lftjz)] ,  (37) 

where  m  is  the  layer  index.  The  characteristic  equation  of  each  layer  Is  the  same 
as  Eq.  (7).  However,  layers  differ  from  each  other  by  different  values  of  r,  which 
is  a  measure  of  the  amount  of  metal  within  a  grating  period.  Following  the  method¬ 
ology  of  the  square-well  grating,  we  match  the  fields  at  the  boundaries  of  each 
region  and  subregion.  Hence,  we  have 


1  ?o(V  *  5 -  5* 

(38) 

r1*n(|n,nt1?„+|n!n+l5n>  '  *!• 

(39) 

fV*o-5o> * 

(40) 

~  0n^5n,n+l-n”  £n,n+l?n^  “ 

(41) 

~m-l  ^~m-l,m\-l +  Im-l,m|m-l  ^  *  ?m+l^|m,nAn+  ?m,m?m^ 

(42) 

9m-l^~m-l,m*m-l  ”  £m-l,m?m-l)  “  ~m+l  ^m.rrAm’  ~m,m?m^  » 

(43) 

with  m  =  1,  ...  n,  and  the  grating  has  n+1  layers  with  a  thickness  of  h.  The 
elements  of  the  l  matrix  are  given  as 

^u.v’p.q  '  Mp('1APVh,{P9-  (44 

All  other  matrices  are  defined  as  In  Part  A  but  with  a  subregion  Index.  It  Is 
trivial  to  generalize  the  formalism  to  treat  a  different  thickness  for  each  layer. 
For  simplicity,  we  choose  each  layer  to  have  the  same  thickness. 


Substituting  Eqs.  (46)  and  (48)  Into  (45),  we  have 


The  expressions  for  R  and  T  are  then  given  as 

R  *  (UK  +  M)  +  UL  +  N)n)}“]{(^K  +  M)  -  (£L  +  N)n}D  (52) 

and 

T  *  A’^L'1  +  N"'*C)”1(L”1K  -  N_1M)(D  -  R).  (53) 

Equations  (52)  and  (53)  are  now  all  we  require  to  calculate  the  fields  at  the  peaks 
and  troughs  of  the  grating. 

III.  Results  and  Discussion 
A.  Square-Well  Grating 

Numerical  calculations  have  been  carried  out  for  a  diffraction  system  with  an 

Incident  wavelength  of  X  *  6471  A,  a  grating  periodicity  of  d  ■  1  pm,  and  the  Incident 

angle  fixed  perpendicular  to  the  grating  surface.  Under  this  configuration,  the 

(first-order)  perturbation  theory  predicts  no  resonant  excitation.  The  diffraction 

amplitudes  are  computed  as  a  function  of  the  groove  depth  for  a  silver  grating  with  the 
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dielectric  constant,  e  *  -17.42  +  0.581,  chosen  to  fit  experimental  data.  Results 
are  plotted  In  Fig.  3  for  the  three  lowest  diffraction  orders.  Since  the 


angle  of  emergence  for  Is  complex,  R2  Is  a  surface  wave.  As  shown  In  Fig.  3,  for 
a  ratio  of  the  grating  amplitude  to  the  period  of  0.283,  the  reflected  beam  Is  at 
Its  minimum,  while  the  first-  and  second-order  diffraction  beams  are  at  their  maxima. 

To  understand  this  behavior,  one  can  draw  an  analogy  to  resonant  scattering 
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theory.  In  particular,  the  bound-continuum  Interaction  problem,  where  a  quasibound 
state  plays  the  role  as  a  coupling  between  the  bound  state  and  the  continuum  states. 

At  the  resonant  energy,  an  Interference  structure,  which  Is  due  to  a  competition 
between  two  equally  possible  pathways,  occurs  In  the  cross  section.  In  the  present 
case,  the  grating  roughness  serves  as  the  coupling  between  the  Incident  wave  and  the 
surface  waves.  At  an  optimum  depth,  the  Incident  wave  couples  strongly  with  the 
surface  waves.  Thus,  the  resonance  occurs.  Moreover,  the  direct  an£  indirect  dif¬ 
fraction  channels  interfere  each  other  and  construct  the  Fano-type  Interference 
structure  in  the  first-order  diffraction  beam.  The  resonance  arises  when  the  cavity 
(the  groove)  has  a  depth  of  approximately  half  of  the  Incident  wavelength,  which  Is 
analogous  to  the  case  of  acoustic  resonances  In  an  open-ended  organ  pipe  occurring 
when  (n+jr)A  equals  the  length  of  the  pipe,  where  n  Is  an  Integer.  In  the  present  case, 
we  believe  that  the  deviation  from  precisely  a  half  is  partially  due  to  the  fact  that 
the  delectric  constant  of  the  Ag  grating  has  a  finite  nonzero  imaginary  part,  such 
that  the  damping  mechanism  leads  to  a  width  as  well  as  a  shift  In  the  resonant 
condition.28  A  calculation  similar  to  the  one  of  Glass,  Weber  and  Mill16  would  be 
helpful  to  further  understand  these  frequency  shifts  and  damping  rates,  which  is 
being  carried  out  In  our  laboratory.  As  far  as  we  know,  this  Fano-type  Interference 
resonance  phenomenon  by  diffraction  from  a  grating  has  never  been  reported  before. 

B.  Levered  Grating 

In  this  part,  we  discuss  the  numerical  applicability  of  the  formalism  derived 
In  Part  B  of  Section  II.  We  test  the  formalism  by  separating  the  square-well 
grating  Into  a  fictitious  layered  grating,  as  shown  In  Fig.  4.  We  then  apply  the 
multilayered  grating  formalism  to  this  fictitious  layered  grating  for  four  and  ten 


layers.  Numerical  results  are  compared  with  those  of  the  square-well  grating  in 
Table  I  for  the  first  three  orders  of  diffraction  beams  at  various  grating  depth. 
First,  the  results  for  four  layers  are  almost  identical  with  those  for  ten  layers. 
This  implies  that  the  method  is  stable  for  a  multilayered  configuration,  which  is 
Important  for  modeling  gratings  of  different  shapes.  Second,  the  layered  grating 
results  are,  in  general,  in  good  agreement  with  the  square-well  grating  results. 


except  near  resonance,  especially  when  the  numbers  are  small.  We  believe  this 


discrepancy  is  due  to  the  degree  of  accuracy  of  the  matrix  inversion,  since  In  the 


layered  grating  formalism  one  needs  to  Invert  matrices  which  are  twice  the  size  of 


those  in  the  square-well  grating  approach.  These  matrices  appear  to  become  singular 


whenever  the  scattering  is  at  resonance.  Moreover,  we  find  that  the  expression  for 


Tn,  Eq.  (53),  is  unstable,  especially  when  the  grating  amplitude  is  large.  This  also 


occurs  in  the  square-well  case  when  Tn  takes  the  form  of  Eq.  (53).  A  more  stable 


expression  might  be  achieved  by  first  evaluating  the  A  and  B  coefficients. 


We  now  turn  to  a  calculation  on  a  three-layered  grating.  The  surface  profile 


has  a  structure  somewhat  between  a  sawtooth  grating  and  a  sinusoidal  grating,  as 


shown  in  Fig.  2c.  We  have  calculated  the  average  field  at  the  peak  and  at  the 


trough  of  the  grating  as  a  function  of  the  grating  depth,  with  the  same  scattering  para¬ 
meters  as  above  and  with  values  of  r  (see  Fig.  1)  as  0.1,  0.3  and  0.5.  The  results  are 


presented  in  Table  II.  We  find  that  Tn  begins  to  exhibit  unstable  behavior  as  the  ratic 


of  the  grating  amplitude  to  the  period  reaches  0.2.  Otherwise,  both  fields  are  con¬ 


vergent  and  stable. 


The  difficulties  in  this  layered  grating  approach  arise  from  the  evaluation 
of  the  nonlinear  eigenvalues  of  Eq.»(7).  As  mentioned  by  Sheng  et  al,i8  eigenfunc¬ 


tions  whose  eigenvalues  form  complex  conjugate  pairs  should  not  be  separated.  While 


truncation  of  the  matrices  is  inevitable  and  the  number  of  eigenfunctions  used  for 


each  layer  should  be  equal,  a  tremendous  amount  of  manual  labor  and  computer  time 


are  required  to  match  these  two  criteria.  This  is  a  potential  problem  in  modelling 


a  realistic  profile. 
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IV.  Summai 


A  new  diffraction  anomaly  is  seen  to  occur  when  the  grating  amplitude  is  approxi¬ 


mately  equal  to  half  of  the  incident  wavelength.  When  the  resonance  condition  is 


met,  the  radiated  energy  (the  scattered  light)  changes  its  "preferred"  direction. 


and  at  the  same  time  surface  waves  are  excited. 


A  new  multilayered  grating  method  has  been  formulated  to  model  grating  pro¬ 


files  of  arbitrary  shape.  The  field  at  the  peak  remains  stable  as  the  grating  depth 
Increases,  while  the  field  at  the  trough  (inside  the  metal)  does  not.  Difficulties 


are  found  in  choosing  an  optimum  basis  set.  A  supercomputer  would  be  helpful  for 


applying  this  multilayered  grating  formalism  to  model  specific  surface  profiles. 
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Table  I 

The  diffraction  amplitudes  of  toe  first  three  orders  for  different  numbers  of 
layers  In  the  square-well  grating  and  for  different  ratios  of  the  amplitude  to 
the  period. 


1 

#  layers 

i 

h 1A 

1 

4 

10 

B 

.10 

.9712 

.9739 

.9739 

1  iR°' 

.20 

.9551 

.9588 

.9588 

.28 

.2367 

.1969 

.1969 

i 

.10 

.1744 

.1567 

.1567 

a  lRi  1 

.20 

.2298 

.2110 

.2110 

.28 

.8551 

.8478 

.8478 

i 

.10 

.0251 

.0310 

.0310 

I  |RJ 

.20 

.0411 

.0430 

.0430 

fl  2 

.28 

.2775 

.2055 

.2054 

y 
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1.  Geometry  and  coordinate  system  for  a  p-wave  Incident  on  a  square-well  grating. 

2.  (a)  Sinusoidal  grating.  The  hatched  area  represents  the  metal. 

(b)  Square-well  grating  showing  a  separation  Into  three  layers,  one  of  which 
Is  periodic  In  the  *- direction  and  two  of  which  are  uniform. 

(c)  Generalization  of  the  square-well  grating  In  which  there  are  three  periodic 
layers. 

3.  (a)  The  magnitudes  of  the  zeroth-  and  second-order  diffraction  amplitudes  plotted 

as  a  function  of  the  ratio  of  the  amplitude  to  the  period  (d  »  1000  ran). 

(b)  The  magnitude  of  the  first-order  amplitude  plotted  as  a  function  of  the  ratio 
of  the  amplitude  to  the  period  (d  *  1000  ran),  with  r  *  0.5  . 

4.  Schematic  diagram  for  a  fictitious  square-well  layered  grating. 
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